Phantom cosmology as a simple model with dynamical complexity 
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We study the Friedmann-Robertson- Walker model with phantom fields modelled in terms of scalar 
fields. We apply the Ziglin theory of integrability and find that the fiat model is non-integrable. 
Then we cannot expect to determine simple analytical solutions of the Einstein equations. We 
' demonstrate that there is only a discrete set of parameters where this model is integrable. For 

, comparison we describe the phantoms fields in terms of the barotropic equation of state. It is shown 

■ that in the contrast to the phantoms modelled as scalar fields, the dynamics is always integrable 
^S) ' and phase portraits are contracted. In this case we find the duaUty relation. 
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'sf ■ I. INTRODUCTION 

> ■ 

I The recently available measurements of luminosity distances of the type la supernova (SNIa) as a function of redshift 
fSj , have shown that current Universe is in an accelerating phase due to unknown form of repulsive energy 0, Q . The 
T— I ' most popular candidate for this dark energy is the cosmological constant. On the other hand the results of large-scale 
structure surveys and results of measurements of masses of galaxies give best fit for density parameter for matter 
^^m,o = 0.3 0,01 (for review of cosmological parameters see Ref. @). Combining data from SNIa with measurements 
of the cosmic microwave background (CMB) radiation, we obtain SIa.q — 0.7 as the best fit value. The sum of the 
densities fitotai.o = 1-02 ± 0.02 obtained by the Wilkinson microwave anisotropy probe (WMAP) agrees with value 
Qh| predicted by inflation and suggests that our Universe is almost fiat on the large scales. Therefore, the assumption of 
Q . flat model with the cosmological constant is in good agreement with observations. 

}—{ The acceleration of the Universe can be explained in two-fold manner. In first approach it is postulated that there 

c/3 , is some unknown exotic matter which violates the strong energy condition p + 3p > 0, where p is the pressure and 
p is the energy density of perfect fluid. This form of matter is called dark energy. In the past few years different 
^ , scalar field models like quintessence and more recently the tachyonic scalar field have been conjectured for modeling 
the dark energy in terms of sub-negative pressure p > —p. A scalar field with super-negative pressure p < —p called a 
phantom field can formally be obtained by switching the sign of the kinetic energy in the Lagrangian for a standard 

■ scalar field. For example in the Friedmann-Robertson- Walker (FRW) model the phantom field minimally coupled to 
a gravity field leads to p -\- p = —0^, where p^ = —l/24>^ -\- V{(j)V p^ = —1/20^ — F((/>), and V{(f)) is the phantom 



o 



■ 

X 



potential. Such a field was called the phantom field by Caldwell |3 who proposed it as a possible explanation of the 
observed acceleration of the current universe when flmfi ^ 0.2. Note that a coupling to gravity in the quintessence 
models was also explored 's'l. 

The second approach called the Cardassian expansion scenario has recently been proposed by Freese and Lewis 
as an alternative to dark energy in order to explain the current accelerated expansion of the universe. In this scenario 
universe is flat and matter dominated but the standard FRW dynamics is modified by the presence of an additional 
term p" such that = pcS = P + ^Bp"^ , where H = {d\na)/dt is the Hubble parameter; and a is the scale factor. 
However, let us note that this additional term can be interpreted as a phantom field modelled by the equation of state 
p = p{p) = \n{\ + 7) — where p = Pm,oCL~^'^^^^^ ■ Therefore for dust matter we obtain p = {n — and n < 
leads to the phantom field. 
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The phantom scalar field can also be motivated from S-brane in string theor y fiol ITlL ll^ . The non-canonical 
kinetic energy also occurs in higher-order theories of gravity and super-grayity |l3lll4|. At first, phantom fields were 
introduced by bulk viscosity effects which can be present in FRW cosmology. They are equivalent to effective pressure 
Poff = P — 3fi/ , where ^ is a bulk viscosity coefficient. It is because dissipation in general relativity is connected 
(in contrast to friction in classical mechanics) with creation of the energy in the expanding universe by the negative 
pressure contribution 0, 0| . 

Without making some specific assumptions on w{z) it very difficult to constrain it from the SNIa data Because 
the astronomical observations do not seem to exclude the phantom fields which violate the weak energy condition, 
it is interes ting to inve stig ate the theoretical pos sibility to describe dark energy in terms of a phantom field 0, 0, 
ES El m Ea il m Ea El m El ES El lH El El El- The Cardasslan expansion with n < which can be 
interpreted as the phantom fiuid effect is also statistically admissible from SNIa data observations |3^. 

In this paper we ask what kind of dynamics can be expected from the FRW model with phantom field. It is 
well known that the standard FRW model reveals some complex dynamics. The detailed studies gave us a deeper 
understanding of dynamical complexity and chaos in cosmological models and resulted in conclusion that complex 
behavior depends on the choice of a time parameterization or a lapse function in general relativity "sT] . Castagnino 
et al. [Ill showed that dynamics of closed FRW models with conformally coupled massive scalar field is not chaotic 
if considered in the cosmological time. They showed that for all initial conditions the universe will collapse in finite 
time and then conclude that there is no chaos in the model. In their work the monotonously growing function is 
defined along a trajectory which diverge at infinity for arbitrary initial conditions. The same model was analyzed in 
the conformal time by Calzetta and Hasi [s^ who presented the existence of chaotic behavior of trajectories in the 
phase space. 

For the cosmological FRW model with a scalar field the kinetic energy form is indefinite, therefore, the domain 
admissible for motion is i?". The similar situation happens in the Bianchi IX in which, as it was proved by Cushman 
and Sniatycki |40j, trajectories have no recurrence property. 

The standard methods of chaos investigation can be also applied to the wide class of a relativistic system 41]. 
Motter and Letelier explained that this contradiction in the results is obtained because the system under consid- 
eration is non-integrable. Therefore we can speak about complex dynamics in terms of non-integrability rather than 
deterministic chaos. The significant feature is that non-integrability is an invariant evidence of dynamical complexity 
in general relativity and cosmology [isL lil . liM li^ . 

There are different motivations to study noningrability in general relativity and cosmology. One reason is the 
possible physical implications of existence of complexity in the systems which for example could help to explain the 
formation of structures. Another reason is to develop suitable tools to study relativistic system. The next motivation 
is to understand the ultimate implication of the time reparameterization. 

For the FRW model with phantom it can be shown that there is a monotonous function along its trajectories 
and it is not possible to obtain the Lyapunov exponents or construct the Poincare sections. Therefore we turn to 
study of non-integrability of the phantom system and set it in a much stronger form by proving that the system 
does not possesses any additional and independent of Hamiltonian first integrals, which are in the form of analytic 
or meromorphic functions. Of course, it is not the evidence of sensitive dependence of solution on a small change of 
initial conditions. However, it is the possible evidence of complexity of dynamical behavior formulated in an invariant 
way. We study non-integrability in the FRW model with phantom fields and find that non-integrability is a generic 
feature of this model and favors rather non-analytical forms of the equation of state. 

It is useful to distinguish between solvability and integrability. While integrability is intrinsic property of the system 
which impose the constraints on the solutions in the phase space, the solvability is related to the existence of closed 
form solutions [4^ . In this paper we concentrate on first integrals rather than solutions of a system. 

We study nonintegrability instead of chaos because this criterion is invariant with respect to time reparameterization. 
Note that while this program of nonintegrability investigation was explicitly formulated by Motter and Letelier [42j |. 
this idea was materialized in papers by Maciejewski and Szydlowski mi3lilE3- Maciejewski and Szydlowski 
also showed that the Bianchi VIII and IX are noningrable in the sense of non-existence of additional analytic first 
integrals and that the Bianchi VIII model is non-intcgrablc in the sense of nonexistence of mereomorphic first 
integrals |49j. The mereomorphic function possesses only poles as its singularities; roughly speaking it is the quotient 
of analytic functions. The latter method is used in this paper. Ziglin proved independently non-integrability of FRW 
closed model with scalar field in the sense of nonexistence of additional mereomorphic first integrals 50]. In turn 
Morales Ruiz and Ramis proved nonintegrability of the Bianchi IX in the same sense . 

For comparison we consider the FRW model with phantom given by the barotropic equation of state which violates 
the weak energy condition. We obtain that this model is integrable in contrast to the previous treatment of phantom 
cosmology. Assuming the barotropic form of the equation of state for the phantom model we obtain the integrable 
dynamics at very beginning. 
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II. HAMILTONIAN DYNAMICS OF PHANTOM COSMOLOGY 



We assume the model with FRW geometry, i.e., the hne element has the form 

ds2 ^ a^(^rj)[-dT]^ + dx^ + fix){d9^ + sin' Od^^)], (1) 

where 

{sinx, < X < fc = +1 
X, 0<X<^k = (2) 
sinhx, 0<x!iOo fc = — 1 

fc = 0, ±1 is the curvature index, < ^5 < 27r and < 6 < n are comoving coordinates, rj stands for the conformal 
time such that dt/a = drj. 

It is also assumed that a source of gravity is the phantom scalar field ip with a generic coupling to gravity. The 
gravitational dynamics is described by the standard Einstein-Hilbert action 

Sg d'x^giR - 2A), (3) 

where — (SttG)^^; for simplicity and without loss of generality we assume 47rG/3 = 1. The action for the matter 
source is 

Sph = -l f d'^x^i-g'^-'ij^.i^, + 2U{^) + iR^^). (4) 



2 . 

Let us note that the formal sign of Ij-iAlP is opposite to that which describes the standard scalar field as a source of 
gravity, where U{tp) is a scalar field potential. We assume 

U{^:) = imV' + ^A^4 (5) 

and that conformal volume J d^x over the spatial 3-hypersurface is a unit. ^ is a coupling constant of scalar field to 
the Ricci scalar 

„ / a fc \ 

If we have the minimally coupled scalar field then ^ = 0. We assume a non-minimal coupling of the scalar field 
. . 

The dynamical equation for phantom cosmology in which the phantom field is modelled by the scalar field with 
an opposite sign of the kinetic term in action can be obtained from the variational principle 8{Sg + 5'ph) — 0. 
After dropping the full derivatives with respect to the conformal time we obtain the dynamical equation for phantom 
cosmology from variation (5(5g+S'ph) /5g = as well as the dynamical equation for field from variation i5(5'g+5ph) / 5^ = 


4^ + 3Hij; = ^ + C-R^- (7) 

dip 

It can be shown that for any value of ^ the phantom behaves like some perfect fluid with the effective energy and 
the pressure in the form which determines the equation of state factor 

Formula ^ differs from its counterpart for the standard scalar field js^ by the presence of a negative sign in front of 
the term tp'^. 

The second derivative {tp'^" in the expression for the pressure in eq. ijHl can be eliminated and then we obtain 

P0 - (-^ - 2^) + ^H{i,'y + 2m - 1)^^' + - - U{i>) + 2^V^- (9) 
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Of course such perfect fluid which mimics the phantom field satisfies the conservation equation 

p^+3H{p^+p^,) = 0. (10) 

We can see that complexity of dynamical equation should manifest by complexity of . 

Let us consider the FRW quintessential dynamics with some effective energy density given in eq. (jSJ. By the 
quintessence we usually understand models with dark energy consisting of a dynamical cosmic scalar field. This 
dynamics can be reduced to the form like of a particle in a one-dimensional potential [s^ and the Hamiltonian of the 
system is 

H{d,a)^ — + V{a) = 0, V{a) ^ -p^^a^. (11) 

The trajectories of the system lie on the zero energy level for flat and vacuum models. Note that if we additionally 
postulate the presence of radiation matter for which p^ cx then it is equivalent to consider the Hamiltonian on 
the level Ti. = E ~ const. Of course the division on kinetic and potential parts has only a conventional character and 
we can always translate the term containing ip"^ into a kinetic term. 

Let us consider now both case of conformally and minimally coupled phantom fields. 



A. Conformally coupled phantom fields 



For conformally coupled phantom fields we put ^ = 1/6 and rescale the field ip (j) = ^a. Then the energy function 
takes the following form for simple mechanical system with a natural Lagrangian function C = l/2gaf)cfq^ ~ ^(9) 

In contrast to the FRW model with conformally coupled scalar field the kinetic energy form is positive definite like for 
classical mechanical systems. The general Hamiltonian which represents the special case of two coupled non-harmonic 
oscillators system is 

n - + V{q) - ^(P' + vD + + By" + Cx^ + Dy^ + Ex'y' , (13) 

where A,B,C,D, and E are constants. 



B. Minimally coupled phantom fields 

For minimally coupled phantom fields (^ = 0) the function of energy takes the form 

£ = ^ + i(0a - c^af - ^a' - - ^-m^c^'a' (14) 

where pcff = -1/2-02 + Uiip), V = -pcsa'^, H = Ijlb? + y(a, V^), = ai\), U{->jj) = 1/2to2^2 i/^x^-i [g assumed. 
This time we parameterize the dynamics by taking variable ip in the original cosmological time and the Lagrangian 
function takes the following form 

a'^ a^d)'"^ lT,,l.,l.r, , , 

2 2 2 4 4 ^ ^ 

where the prime denotes the differentiation with respect to the cosmological time parameter t, and V = —posa^, 

Peff=-l/20'2-|-C/(0). 



III. NON-INTEGRABILITY AS AN INVARIANT FEATURE OF PHANTOM COSMOLOGY WITH 

SCALAR FIELD 

For a given Hamiltonian system, it is difficult to show that the system under consideration is non-integrable. In 
general, there are two formulations of necessary conditions for the integrability presented by Ziglin |54l l55j | and 
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Morales- Ruiz and Ramis |56l l57j . Both approaches base on a deep connection between properties of solutions in an 
enlarged complex time plane and the existence of first integrals. This idea originates from works of Kovalevskaya and 
Lyapunov. 

Let us study the general case of the Hamiltonian system describing the conformally coupled phantom field in the 
FRW model of the universe. We have 

n^l{pl+pl) + Viq,,q2), (16) 



2 

V{qi,q2) = ^ 



^ 4 I 4 2 2 2 



where qi — a, q2 = 4>, pi — a, p2 = (j), A = —A and A = —A. This Hamiltonian has the natural form in which the 
potential is a homogeneous function of degree four with respect to both variables a, 4>. 

From the point of view of complex dynamical behavior it is useful to distinguish from Hamiltonian (|16|l with 
w? = —y? < 0. This case is interesting because of spontaneous symmetry breaking |58]. The Poincare section in this 
case can be obtained as well as the Lyapunov exponents. In this model the chaotic behavior is present. 

In other cases we can define by analogy to Castagnino et al. the monotonic function along trajectories. From 
this fact we obtain that trajectories escape to infinity and the system has no recurrence property which guarantee the 
topological transitivity (the standard chaos indicators cannot be obtained). 

Motter and Letelier ar gue d that the cosmological systems with scalar fields are non-chaotic but complex in the 
sense of non-integrability j42l | . Moreover the non-integrability is an invariant property of system under the coordinate 
change. 

In the second distinguished case the complexity has the same character, and we apply the some tools to confirm 
the Liouville non-integrability of this system. The Liouville integrability of the Hamiltonian system means that there 
is as many functionally independent functions which are in involution (Poisson brackets vanish) as is the dimension 
of the system. 

Now we consider the problem of non-integrability in both cases. The non-integrability of the non-flat first case with 
the spontaneous symmetry breaking was investigated by Ziglin ( A = A = — the Yang-Mills potential) ''so] . In turn, 
we apply the Ziglin and Morales-Ruiz and Ramis methods to flat phantom models with conformally coupled scalar 
flelds with arbitrary parameters. 

The integrability of Hamiltonian systems with a natural Lagrangian was analyzed in details by Yoshida Isol Isol Isil 

t ls^ in the framework of Ziglin's approach. Later Yoshida's results were sharpened by Morales-Ruiz and Ramis 
. Note that we applied the Morales-Ruiz and Ramis result to system H16(l . but with the indeflnite kinetic energy 
formT= l/2((pf-pi). 

The counterpart of Hamiltonian ()16|l for a standard scalar flcld can be obtained after the canonical transformation 
of variables 



and 



qi^Qi, pi^Pi: 



92 iQ2, P2 ^ P2 = -iP2- 



Then of course dp2 A dq2 = dP2 A dQ2. However, in this case the phase space is complex. Moreover, trajectories 
have no recurrence property which guarantee the topological transitivity, which an essential element of the standard 
understanding of chaos. 

The fundamental papers of Ziglin |54l Issj gave the formulation of a very basic theorem about non-integrability 
of analytic Hamiltonian systems. The Ziglin idea connects properties of solutions on a complex time plane and the 
existence of first integrals. This approach takes its origins in works of Kovalevskaya and Lyapunov. 

The Yoshida criterion is presented in Appendix. We apply this criterion to the analyzed system. Then the equation 

q^V'{q), 9= (91,92) 

has the following solutions 



zi 



The integrability indices for this points are 

A, = -try"(z,) - 3, i = 1,2,3 
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and 




AiA2-2(Ai + A2) + 3 
I-A1A2 




(17) 



Thus, from the Yoshida criterion foUows that if there exists / e {1,2,3} such that A; € A'4 then system (I16f) has no 
additional meromorphic first integral that is functionally independent of H . Moreover, our previous application of the 
Morales-Ruiz and Ramis result to the considered system gives that if we introduce quantities H17() and three discrete 
sets 



then if Ai , A2 , A3 ^ / = Ii U /2 U /a the system is non-integrable. Therefore, only for certain values of model parameters 
the phantom cosmology is integrable. We can conclude that the Liouville non-integrability is the generic property of 
the system. 

If we consider a non-flat model then the effects of curvature are negligible near the singularity and the considered 
case describes a generic situation. In this way the phantoms give rise to the complex dynamics in the sense of non- 
existence of a sufficient number of independent first integrals. As a consequence, we can express some scepticism 
about prediction for the equation of state factor w{z) in the presence of the phantom component of dark energy. 

Our conclusion is that in a generic case the phantom scalar field can produce the complex behavior. The complexity 
of dynamics is formulated in terms of non-integrability (i.e., non-existence of an additional first integral) because a 
standard understanding of chaos has no significant physical meaning in the context of a gauge freedom in the choice 
of a lapse function (time parameterization). 

Beck proposed an interesting idea that stochastically quantized scalar fields can offer some solution to the 
cosmological coincidence problem of A. In this approach the chaotic fields have a classical equation of state close to 
p = — p, i.e., that the chaotic fields naturally generate a small cosmological constant. It is possible that phantoms are 
just a phenomenological description of this situation on a purely classical level. 

Let us also note that regular behavior of dynamics in phantom cosmology can appear to be different from the 
considered types of potentials (6^ but ¥{({)) oc (f^ is the simplest one in which this phenomenon occurs. Moreover it 
can only appear if we treat phantom energy in terms of a single scalar field. 

IV. PHANTOM COSMOLOGY IN TERMS OF BAROTROPIC EQUATION OF STATE VIOLATING 

THE WEAK ENERGY CONDITION 

It is well-known that for given evolution of the model it is possible to construct a potential for a minimally coupled 
scalar field which would reproduce this cosmological evolution 66] . Sometimes it is possible to find the explicit form 
of scalar field potential can reproduce the evolution arising in some perfect fluid cosmological model '67] . 

The very different picture is found, if we consider phantom energy as a some kind of perfect fluid with super-negative 
pressure then, in the contrast to previous case, the dynamics is regular at very beginning. 

Let us consider the dynamics of the FRW models with phantoms where the specific form of the equation of state 
for phantom fluid is assumed. We model the fluid which violates the weak energy condition using the equation of 
state p = wp and w — const < — 1. Such a model of fluid can be treated as the simplest phenomenological model of 
phantom matter. 

The dynamics of this model can be represented by a two-dimensional dynamical system (therefore non-chaotic at 
very beginning) on the phase plane {x, x) = {x, y) or by motion of a classical particle in the one-dimensional potential 



h = {p{2p-i) \pez}, 

/2 = {1/8 [-1 + 16(1/3 + j5)2] |peZ}, 
/3-{l/2[3/4 + 4p(p-l)] |peZ}, 



(18) 



V{x): X = a/ao H^l, i-e. 



X 



y, 



y 



dV 
dx ' 



(19) 



System H19|) has the first integral of energy in the form 




(20) 
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where 

V{x) = (an,o.T-i + r!ph,o2;~3(i+"')+' + r!fc,o) ■ (21) 

For the mixture of noninteracting matter and phantoms here O^^o are the density parameters at the present epoch. 
In the general case the potential of the particle-universe takes the form 

y(2;)^-i^^],ox-=^(l+-)+^ (22) 

i 

where Wi — —1 for the cosmological constant, Wi = —1/3 for string fluid (also curvature fluid), Wi — —2/3 for 
topological defects. 

The phase portraits for the model described by system H19|) for the potential function (|21|) are shown on Fig. ^ 
The trajectory of the flat model separates the regions of closed (Ofc,o < 0) and open {Q.kfi > 0) models. Moreover, 
both phase portraits are topologically equivalent. The presence of additional terms like strings, topological defects 
(see [31|) do not change the structure of the phase plane. There is the single critical point located on z-axis as an 
intersection with the boundary of the strong energy condition p + 3p > 0. 

Note that the obtained phase portraits are equivalent to phase portraits of the FRW model with the cosmological 
constant. 




FIG. 1: The phase portraits for the FRW models with phantom matter described by the equation of state p — wp for (a) 
w = —4/3 and (b) w — —5/3. The dashed lines are the flat model trajectories. The shaded region is the region of accelerated 
expansion of the universe. Note the topological equivalence of both phase portraits. 

From the first integral H20|l for the mixture of the cosmological constant and phantom type matter in the flat FRW 
models we obtain the relation 

(inx) ' = (^) = r!ph,ox-3(i+») + Oa,o, (23) 

which preserves its form structure under the change both a position variable and a sign of the quintessential parameter 
{w + 1) 

x^-, (l + w) ^ -(l + w). (24) 

X 

Therefore, if a phantom epoch exists its dynamics can replicate the corresponding evolution for the sub-negative 
equation of state (for example w — —4/3 corresponds to w ^ —2/3). From this kind of symmetry we obtain that if 
x{t) is the solution of l|23() for the sub-negative equation of state p — wp then x~^{t) is also its solution for other form 
of the negative equation of state p = —{w + 2)p. Let us note that for w = —1 the duality relation which is motivated 
by superstring theory of duality symmetries |68l l69j| is the exact symmetry of dynamical equations. 
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V. CONCLUSIONS 



This paper addressed the problem of complexity of the flat FRW dynamics with phantom modelled in terms of scalar 
fields. We proposed a criterion of non-integrability in the Liouville sense as an adequate measure of complexity of the 
phantom cosmology. This approach is opted because the gauge freedom in the choice of a time parameterization or a 
lapse function unables us to discuss chaos in general relativity in the standard way. Non-integrability is an invariant 
feature of a system and can be use as an indicator of its complex dynamics. 

We considered the two approaches to model the phantom fields in the FRW model and showed how different 
dynamics of the models are in these approaches. In the first approach which could be called "microscopic" we used 
scalar fields in modelling the pressure and energy density. In the second approach wc used the explicit dependence of 
pressure versus energy density, which can be called phenomenological and found integrable dynamics in contrast to 
the non-integrable dynamics of the first approach. 

In this context non-integrability is a generic feature of the phantom cosmology and only for a certain discrete 
set of values of model parameters phantom cosmology is integrable in the Liouville sense. From the physical point 
of view one can interpret this property as the complexity of dynamical behavior of trajectories in the phase space. 
Therefore we cannot expect any simple analytical relation for solutions (trajectories) of the system or the form of 
relation p{(l), a, (j), a) and a, (j), a) along the trajectories. 

Nevertheless, the integrable cases with zero-measure in the space of all solutions exists, they can lead to some 
analytical dependence of after the elimination of time. As the example of simple analytical form of the equation 
of state we considered the barotropic equation of state which violates the weak energy condition. We obtained that 
this model is integrable and exhibits the regular dynamics. 

We conclude that non-integrability which is a generic feature of the FRW model with phantom fields favors rather 
non-analytical forms of the equation of state. However, assuming the barotropic form of equation of state for the 
phantom model we obtain the integrable dynamics at very beginning. We expect that phenomenological equation 
of state can be realized by a microscopic scalar field with some potential. Therefore, we assume at very beginning 
the existence of some relation (a first integral) between scalar fields, its derivative and evolutional parameter of the 
universe. From the mathematical point of view this requirement means the existence of some invariant in the phase 
space. If we prove non-integrability, than there no such relation in the considered class of potentials. If there is even 
a discrete set of parameters for which the system is integrable we hope for finding this relation. 



APPENDIX: THE OUTLINE OF NON-INTEGRABILITY CRITERION 



Here we present only the facts needed for a formulation of a criterion in possibly simple settings. We consider 
a complex symplectic manifold C^" with the canonical symplectic structure Q,. A Hamiltonian vector field vh is 
determined by a complex Hamiltonian function H : C^" — > C by the relation ^{vht ) = dH . We assume that 
Hamilton's analytic equations 

dz 

-f=VH{z), z = (zi,...,z2„)eC2", tec, (A.l) 
at 

have the non-equilibrium solution z — Lp{t). To simplify the exposition we assume that this solution lies on a two- 
dimensional invariant plane 

H = {(zi, . . . , Z2n) e C^" I = 0, z = 1, . . . , 2(n - 1)}. 

The phase curve F = {'p{t) S C^" 1 1 G C} is a Riemannian surface with a local coordinate t. Together with equations 
IjA.ljl we consider also variational equations along solution ip{t) 

§=Am, A{t) = ^Mt)). (A.2) 

This system separates into the normal and the tangential subsystems. In our settings this separation takes a very 
simple form — the matrix A{t) has a block diagonal structure. We consider the normal variational equations (NVE) 

^^B{t)r^, rjeC'^-~'\ (A.3) 

where B{t) is 2{n — 1) x 2(n — 1) upper diagonal block of matrix A{t). We choose a point G C and a matrix of 
fundamental solutions of the NVE X{t), defined in a neighborhood of to. With a close path a on complex time plane 
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starting and ending at point to we can associate a matrix S G GL(2(n — 1),C) in the following way. We integrate 
NVE lfO|) along the path a, i.e., we make an analytic continuation of X{t) along this path. As a result from the 
fundamental solution X{t) we obtain another fundamental solution Y{t). From the general theory of linear systems 
it follows that Y{t) = SX{t) for some S € GL(2(n — 1),C). Because the system is Hamiltonian, S' is a symplectic 
matrix, i.e, S € Sp(2(r7, — 1),C). In this way, considering all possible paths, we obtain a matrix representation of 
the first homotopy group 7ri(r) of T. It forms a finitely generated subgroup of Sp(2(n — 1),C) and it is called a 
monodromy group. We denote it M . 

Let us take an element of the monodromy group g e M . Its spectrum has the form 

spectr(g) = (Ai,Aj;\...,A„_i,A;^:^i), Ai e C. 

The element g is resonant if 

n-l 

Y[ Af' = 1 for some (fci, . . . , kn-i) G Z"-i\{0}. 
1=1 

Theorem 1 (Ziglin (H3|) Let us assume that there exists a non-resonant element g G M . If the Hamiltonian system 
possesses in a connected neighborhood ofT n — 1 meromorphic first integrals which are functionally independent with 
H then for an element g' G M: if ge — Ae for A G C and e G C^^"~^\ then g{g'e) = X'(g'e) for some A' G C. 

In the case of a system with two degrees of freedom this theorem can be formulated in a more operational way. 

Theorem 2 Let us assume that there exists a non-resonant element g G M . If there exists other element g' G M 
such that 

1. trg' 7^ and gg' ^ g' g, or 

2. trg' = and gg' g ^ g' , 

then there is no additional meromorphic first integral functionally independent of H in a connected neighborhood of 

r. 

The main difficulty with the application of the Ziglin theorem is the determination of the monodromy group of the 
NVE. Only in very special cases we can do this analytically. Yoshida |^ developed the Ziglin approach 

for these cases when the Hamiltonian of a system has the natural form and the potential is a homogeneous function. 
In this case a particular solution can be found in the form of 'straight line solution' and the NVEs for it can be 
transformed to a product of certain copies of hyper-geometric equations for which the monodromy group is known. 
This allows to formulate adequate theorems in a form of an algorithm. Below we describe it for the Hamiltonian 
system with two degrees of freedom. 

Consider the Hamiltonian 

H=l{pl+Pl) + V{qi,q2), ((?l,(?2,Pl,P2)GC^ (A.4) 

where V{qi,q2) is the homogeneous function of degree fc, i.e., 

V{Cqi,Cq2) = C''V{quq2). (A.5) 
In a generic case this system has a straight line solutions of the form 

gi = C'Mt), 92 = C20(i) (A.6) 
where (j){t) is a solution of a nonlinear equation 

and (Ci, C2) 7^ (0,0) are solutions of the following system 

Ci = diV{Ci, C2), C2 = d2V{C\, C2). (A.7) 
The variational equations take the form 



. ... 




>ii 


V12' 










V21 


V22 
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where Vij = didjV{Ci,C2) for i,j = 1,2. Since the Hessian of V is symmetric it is diagonaHzable by an orthogonal 
transformation and the system separates to 



(A.8) 
(A.9) 



where Ai, A2 are real eigenvalues of the Hessian. Let us note that it is not true for indefinite systems where the 
Hessian is not a symmetric matrix. 

It can be shown that the Hessian of y at C = (Ci, C2) has the eigenvalue Ai — k — 1. Thus, its second eigenvalue 
is equal A2 — try(Ci,C2) — (fc — 1), and it is called the integrability index. Equation (|A.9|) can be transformed to 
the hyper-geometric equation. The monodromy matrices of this equation are parameterized by A and conditions of 
the Ziglin theorem put a restriction on values of A — simply, we can identify those values of A for which the system is 
not integrable (more precisely: it does not possess an additional meromorphic first integral). To state it accurately 
let us define 



hip) 



\p{p + 1) 



k 

P, -^P{P 



I) 



p+l 



peN, 



(A.IO) 



and 



iVfc =R \ U hip). 

peN 



(A.ll) 



Then it follows that Hamiltonian system (jA.4p with homogeneous potential ljA.5|) of degree k is not integrable if the 
integrability index A corresponding to a certain straight line solution ljA.6|) belongs to Nk- Let us note that equations 
(|A.7p usually have several solutions and thus it is necessary to check the Yoshida criterion for each of them. 
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